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Abstract 

£h A general recipe is developed for the study of rigid body dynamics 

in terms of Poincare surfaces of section. A section condition is cho- 
sen which captures every trajectory on a given energy surface. The 
possible topological types of the corresponding surfaces of section are 
t-H determined, and their 1:1 projection to a conveniently defined torus is 

proposed for graphical rendering. 
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1 Introduction 

O 

The dynamics of a rigid body in a constant field of gravity, with one point 
J^j held fixed, is in general non-integrable. The configuration space S0(3) has 

three degrees of freedom, but except for the special cases of Euler, Lagrange, 
and Kovalevskaya, the system possesses only two integrals of motion: the 
energy H = h and the angular momentum component in the direction of 
gravity, L z = I. Liouville integrability would require a third integral, but 
this does not exists. The parameter space is essentially 4-dimensional: when 
lengths, times, and energies are properly scaled, there remain two freedoms 
for the principal moments of inertia A = diag( J 4i, Ai, A3), and two for the 
location of the body's center of mass r = (r\,r2,rz) relative to the fixed 
point. Within this 4-D set of parameters, the Euler case r = defines 
a 2-D subset (two freedoms in the moments of inertia), Lagrange's case a 
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1-D subset (one free ratio of the moments of inertia, no freedom for the 
direction of r), and Kovalevskaya's case is a single point (no freedom in the 
moments of inertia nor in the direction of r, provided the body-fixed frame 
of reference is properly chosen). This implies that the vast majority of rigid 
body systems exhibits some degree of chaotic motion which to this very day 
has largely remained unexplored. 

The three degrees of freedom of configuration space SO(3) may be re- 
duced to effectively two in view of the S ^symmetry with respect to the axis 
of gravity. Thereby the configuration space reduces to the Poisson sphere 
S 2 , and the corresponding 4-D phase spaces T^S 2 are labeled with the an- 
gular momentum constant /. The surfaces £?, of constant energy H = h are 
manifolds in T^S 2 , except at values (h,l) where their topological character 
undergoes a bifurcation. The first step in the analysis of phase space struc- 
ture is the identification of the topology of and of its bifurcation scheme. 
This requires to study the energy-momentum map from phase space to the 
(h, Z)-plane, and to determine its critical values. Already this first step is 
more difficult than one might think. It has been worked out for certain sub- 
sets of parameters [7, 10] but not for the entire 4-D family of rigid bodies. 
However, even though the bifurcation schemes have not been resolved for all 
cases, it is known that the connected components of S\ , come only in four 
types [11]: sphere S 3 , direct product S 1 x S 2 , real projective space MP 3 , or 
connected sum (S 1 x S 2 )#(S 1 x S 2 ). 

Once the energy surfaces £\ l are given, the problem is to find how they 
are partitioned into regular and irregular types of motion, i. e., into invariant 
subsets of one, two or three dimensions (isolated periodic orbits, resonant 
or non-resonant tori, and chaotic regions, respectively). The most valuable 
tool for this kind of studies is the method of Poincare sections. It requires 
the identification of a 2-D surface of section V\ i C £ ? , which intersects all 
possible trajectories. Finding such a global Poincare section is a non-trivial 
matter. It is in general not possible to choose a surface which intersects all 
trajectories transversally [2]. However, it is possible to find a surface (or 
a set of disjoint surfaces) which is complete in the sense that every orbit 
intersects - or at least touches - it repeatedly. A constructive procedure to 
obtain such section conditions was given in [5] and will be employed here. 

The next question refers to the topology of the surfaces of section, 
and how the various types of "P 2 1 may be adequately represented in two- 
dimensional plots. This is the main concern of the present paper. We show 
that a single connected component of V\ { may be either a sphere S 2 , a 
torus T 2 , or a 2-D manifold M 2 of genus g = 2, 3, or 4. A single kind of 
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graphical representation applies to all cases: the 1:1 projection to a torus 
constructed from two copies of a two-fold punctuated Poisson sphere. We 
call this the "PP-torus" T|(7) and propose it as a convenient universal 
tool for investigations of the complex dynamics of rigid bodies. For exam- 
ple, it may be interesting to vary (h, I) for fixed parameters (Ai, A2, A3) 
and r2, rs), in order to obtain a complete survey of the phase space 
structure for a given rigid body. For the integrable Kovalevskaya top such 
a survey was presented in [9] , but the new tool can be applied to any choice 
of parameters. Alternatively, one might want to vary parameters at fixed 
(h,l), and to follow the fate of certain conspicuous features (like isolated 
periodic orbits, or major chaotic regions). 

The paper is organized as follows. Section 2 recalls how the topology 
of energy surfaces is determined from the effective potential on the Poisson 
sphere. The new Poincare section is introduced in Section 3 and compared 
to a proposal made earlier in [6] . Finally, Section 4 describes the topology of 
the surfaces of section V\ t and their projection first to the Poisson sphere, 
then to the PP-torus constructed from it. 

2 Topology of energy surfaces 

The heavy rigid body on T*SO(3) is symmetric with respect to rotation 
about the axis of gravity. Reduction by this symmetry gives the Euler- 
Poisson equations for the motion of a rigid body about a fixed point as seen 
in a co-moving frame, 

7 = {-y,H} = 7 x cj 
i = {1, H} = 1 X u> — mg-f x r . 

Here 7 is the unit vector along the spatial z-axis (the axis of gravity), I 
is the angular momentum vector, u = A~ l l the angular velocity vector, 
A = diag(Ai, A2, A3) is the matrix of principal moments of inertia, and r is 
the position of the center of mass in the body. The connection to Euler's an- 
gles with respect to the z-axis is (71,72,73) = (sin V>sini9, cos ip sin??, cos 1?); 
the angle ip of rotation about the z-axis is eliminated by the symmetry re- 
duction. The remaining angles (1?, ip) parameterize the Poisson sphere S 2 . 
The constants mg can be absorbed into r and will be ignored from now on. 
These equations on the reduced phase space are Hamiltonian with respect 
to the Poisson bracket 

{F, G} = (V-y-F, 7 x ViG) + (VjF, 7 x V^G + 1 x ViG) (2) 
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and Hamiltonian 

H( 7 ,l)= 1 -(l,A-H)-{ 1 ,r), (3) 

where (•,•) denotes the standard Euclidean scalar product in M 3 . The 
bracket has two Casimirs 1 , 

7( 7 ) = <7,7> = 1 and L z (j, I) = (7, 1) = I. (4) 

Fixing the Casimirs to the values = 1 and L z (~f,l) = I defines the 

reduced phase space T^S 2 . Fixing, in addition, the energy H = h defines 
the reduced energy surface S\ v The effective (or amended) potential on S 2 
is, see e. g. [9], 

The term proportional to I 2 (the kinetic energy of rotation about the axis 
of gravity) contains the moment of inertia (7,717). The accessible region 
on S 2 for fixed energy h and angular momentum I is 

U h ,i = {7 : 17,(7) < h) (6) 

The energy- Casimir map is given by 

F : M 6 — >R 3 (h,l,l) 
(7, 1) -^((#(7, 0.^(7, 0,7(7)) 

Its critical points (7, 1) : rankDi 7 < 3, are the relative equilibria 2 , and 
their images under F are the critical values. Since the last component of 
the energy- Casimir map is constant, it suffices to consider the (h,l) energy- 
momentum plane as the image. The energy surface is the preimage of 
the point (h, I). The set of critical values is called the bifurcation diagram. 
For critical values the energy surface is not a smooth manifold, and the 
topology in general changes upon crossing critical values. 

Instead of computing the rank of F, the equilibria of (1) can be computed 
directly. The relation between the two approaches is that the gradients of 
the Casimirs are in the kernel of the Poisson structure, and if rank DF is 
lower than 3, the gradient of the Hamiltonian is zero or a linear combination 
of the gradients of the Casimirs. A third method to obtain the relative 
equilibria is to compute the critical points of the effective potential XJ\. The 



(7) 



1 note that I does not denote the length of I, but the value of its z-component (7,/). 
2 see Arnold [1], Appendix 5C. 
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Figure 1: The bifurcation diagram showing critical values (h,l) of the energy- 
Casimir map F for the Kovalevskaya case A = diag(2, 2, 1), r = (1,0,0). Right: 
Magnification of the rectangle shown on the left. 

energy surface can be viewed as a singular circle bundle over the accessible 
region Uh,i, see e. g. [2]. At critical points of Ui the topology oiUh,i changes, 
inducing a change in the topology of the energy surface. Thus a bifurcation 
diagram as shown in Fig. 1 can either be read as a statement about the 
existence of relative equilibria at the critical values shown, or as a statement 
about the topology of the reduced energy surface for the non-critical values. 

For example [9], along the line I = 1.84 in Fig. 1, the topologies of the 
accessible regions U^i are, from left to right, empty set, disk, disk with 
two holes, disk with one hole, and all of S 2 . The transitions take place 
when crossing the lines of critical values in the bifurcation diagram. The 
corresponding topologies of the energy surface l are determined by the 
topology of the accessible region, see [2]. They are 0, S 3 , (S 1 xS 2 ) # (S 1 xS 2 ), 
S^S 2 , MP 3 , again from left to right. More general cases of the bifurcation 
diagram have been discussed in [7, 11], see also [6]. 

It should be remembered that of the three Euler angles, only #(t) and 
ip{t) appear in the Euler-Poisson equations; the angle ip(t) must be deter- 
mined by integration of [9] 

Hence, relative equilibria 7 = const, I = const of the symmetry reduced 
system correspond to periodic motion of the full system. 
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3 Poincare surface of section 



For systems with two degrees of freedom and 3-dimensional compact energy 
surfaces, the relevant information about phase space structure (stability 
of periodic orbits, relative extent and entanglement of regular and chaotic 
motion) is contained in the 2-dimensional maps induced by the motion on 
suitably chosen Poincare surfaces of section. This has become a standard 
tool from celestial mechanics to molecular dynamics, yet in the study of rigid 
bodies there have so far been only few applications [3, 4, 8]. However, as 
the Euler-Poisson equations describe a system with effectively two degrees 
of freedom (at given I), Poincare sections are the ideal method to study the 
complexity of rigid body dynamics. In the following, we propose a variant 
which is at the same time general and easy to implement. We begin with a 
few general considerations. 



3.1 General features of Poincare sections 

It would be nice to identify a two-dimensional submanifold V\ i of the en- 
ergy surface £f ll in such a way that every orbit meets it repeatedly and 
transversally. As was shown in [5] , complete sections can indeed be defined, 
i. e. sections which capture every single trajectory. However, globally trans- 
verse sections do not exist in general, for topological reasons [2]. Namely, 
when the surface of section V\ l divides E\ l into an "inner" and an "outer" 
part, trajectories that "come in" have to "go out" again, so that part of 
V\ i is traversed by ingoing, another part by outgoing orbits. The boundary 
between the two is then a subset of P?> where the trajectories are tangent. 

A complete surface of section V\ l can be obtained with a recipe from [5] : 
given any smooth function W which maps the energy surface into a bounded 
set, define S := W and take £(7, 1) = as section condition: 

Kf-= {(7,0 5(7,0 = 0} • (9) 

V\ l is a smooth 2-dimensional manifold except when the map 



P^j—R^M, 1,0) 
(7,0 ^ (P( 7 , 0,^(7, 1), Hi), 5(7,0) 



(10) 



has a critical point indicated by rankDP(if, L Z ,I, S) < 4. Clearly the rank 
of P drops when the energy surface is not a smooth manifold, but there may 
be additional singular values. 
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Figure 2: Schematic presentation of the Poincare surface of section V\ l (circle) 
with four trajectories two of which are tangent in the thick dots, and its projections 
onto two copies of the sphere S 2 (7) (vertical bars). Left: defined by the sign 
of S (projections not invertiblc). Right: defined by the sign of A. 

The tangent set, if it exists, is given by those points on V?i for which 
S = 0. Consider for example the Hamiltonian H = \p\ + \p\ + V(q\,q2) 
and the section condition S = q2 = 0. 3 The surface of section is then given 
by V\ = {(qi,pi,P2) '■ 2P1 + \v\ + ^(?1j0) = h}. Its usual representation 
is in terms of its projection to the (gi,pi)-plane which is a 2:l-map. To 
make it unique, only the part qi > is considered; the part qi < is 
ignored because the two parts are related by time reversal and contain no 
independent information. The line of tangency 5 = qi = P2 = lies in the 
projection plane and is in fact the boundary of the energetically accessible 
region {(gi,pi) : \p\ + V(qi,0) < h}. 

The situation is similar but somewhat more complicated in the case of 
rigid body dynamics. Assume we have a section condition 5(7, 1) = 0. The 
corresponding surface V\ t C E\ l "lives" in the 6-dimensional (7, Z)-space, 
7 G S 2 and 1 £ R 3 . A convenient projection, as we shall see, maps it 
onto the Poisson sphere, 7r : V\i — * S 2 (7). Here it is not generally true 
that the boundary of the projection coincides with the points of tangency. 
The situation is rather as shown schematically in the left part of Fig. 2. 
There Vi , corresponds to the circle, and the bars to its right are meant to 
represent two copies of S 2 (7). The four trajectories illustrate where "P 2 ; 
intersects incoming and outgoing orbits (S < and S > 0, respectively). 
In this sketch, the line of tangency reduces to two points which may lie 
anywhere on the circle. Then, if S + (7) were chosen to carry the projection 

3 This section condition is not of the form S = W, hence it need not be complete; there 
may be orbits which never cross the line 92 = in the (qi, Q2)-plane. 
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of the part where S > 0, and likewise S (7) for S < 0, these maps would 
obviously not be 1:1. In this setting, 7 could not serve as a local coordinate 

On the other hand, consider the right part of Fig. 2. There the circle 
representing V\ x is divided by the two solid squares in such a way that the 
two projections to S ± (7) are indeed 1:1. A drawback of this welcome fea- 
ture is that both projections contain incoming and outgoing intersections; 
the points of tangency no longer form the boundaries of the projections. But 
this turns out not to be serious. The important point is that the projections 
S ± (7) provide local coordinates on the two halves oiV\ l . But how is this 
schematic picture to be implemented? It must be possible to uniquely de- 
termine the momenta I on V\ l from the coordinates 7. With V\ l defined 
by H = h, L z = I, and S = 0, the implicit function theorem guarantees that 
this can be done unless 

A(7,0:=det^(fr,L*,S) = 0. (11) 

The relevant division of V\ l into two parts is therefore given by the sets 

V h,l C V Kl : A > and V h,l c V h,l : A < which project to S+(7) and 
S~(7) respectively. The condition (11) contains the definition of the dividing 
line on V^ t . It will be seen in (15) that there may exist other lines in "P^ ; 
where A = 0, but these project only to points in S 2 (7). 

Note that even when the system is symmetric under time reversal (which 
is not the case here, at fixed I 7^ 0), it may not be sufficient to consider 
only one of the two projections. As can be seen in Fig. 3, orbits with 
intersection points close to tangency tend to have incoming and outgoing 
intersections with the same sign of A; hence, in order to capture all orbits, 
both projections are needed. 

3.2 The section condition 

The first general section condition of the type S = W was proposed in [6] 
with W = (1,1). It was motivated by a study of the projections of the 
energy surfaces E\ x to i-space. There the preimages of the envelopes of 
these projections could be characterized by S = 0, or (£,7 x r) = 0, i.e. in 
terms of local extrema of the values of the total angular momentum. Clearly 
this section condition is the adequate choice when the envelopes in Z-space 
are used to represent the Poincare surfaces of section; a number of beautiful 
examples have been worked out in this setting (Gashenenko, private commu- 
nication). However, the complex topology of the corresponding surfaces V\ l 
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Figure 3: Only both copies of S 2 (7) allow a complete characterization of the 
phase space structure since intersection points can be projected to one copy only. 

(being manifolds of genus up to five and even higher) makes it hard to study 
them. Furthermore, the /-envelopes are different for each point (h,l), and 
difficult to parameterize, hence they do not lend themselves for extended 
investigations. An obvious way out of this last difficulty is to project the 
surfaces V\ x , defined by W = (1,1), to the Poisson sphere S 2 (7). However, 
it turns out [6] that this produces artifacts in the projections which have no 
relevance in V\i itself: the condition W = (1,1) is well adapted for Z-space, 
not for 7-space. 

Therefore we propose another choice: 

dW 

W = ( 7 ,r) =}► S(~f,l)= =(A- 1 l,rx 1 ) = 0. (12) 

dt 

This is motivated by the behavior of the Lagrange top where r lies on the 
body axis so that (7, r) oc cost?. The angle $ oscillates between i? m i n and 
$max which makes 1) = a natural section condition. In addition, when 
successive points ($ n , ip n ) (n = 0, 1, 2, ...) of the corresponding Poincare map 
are plotted on S 2 (7) , a natural winding number may be read off in terms of 
the increments of if). The section condition (12) is an obvious generalization; 
it picks out extrema of the projection of the center of mass r to the vertical 
direction. As we will show, the resulting surfaces of section V\ , have genus 
only up to 4, and they have fewer bifurcations than those studied in [6]. 
Moreover, the projections ^('P? 7) onto the Poisson sphere coincide with the 
accessible region Uhi, and the condition A(~f,l) = defines a line which 
projects to the boundary dUh,i- Finally, this projection is simple because 
generically it has exactly two preimages; the exceptions are 7 G dU^i which 
have only one preimage, and 7 || r where the preimage is a circle. 

To prove these statements, we simultaneously solve the equations H = h, 
L = I (which together define the energy surface £| >) and S = 0; i.e., given 



9 



a point 7 on the Poisson sphere, we determine the set of corresponding 
values I such that (7,/) £ V^ t - The explicit calculation is given in the 
Appendix. There we excluded the possibility that 7 and r are collinear, 
r x 7 = 0, so let us start here with this case: 7 = ±r/r =: r. The 
condition S = is then identically fulfilled and gives no restriction on I. 
The condition L z = = I defines a plane in i-space, and H = h the 

ellipsoid (I, A~H) = 2(h + (r, 7)) = 2(h ± r). Three cases are possible: 

(i) The intersection of plane and ellipsoid is a topological circle. 

(ii) They do not intersect at all. 

(iii) The plane is tangent to the ellipsoid in a point I*. 

In the latter case the normal to the ellipsoid must be collinear with 7, i. e. 
A~ l l* = £7 or I* = £Aj. From = I we obtain £ = l/(-f,A-f), hence 

the energy equation becomes 

2{h±r) ^ h = U l { 1 ) = U l {±v). (13) 



(7,^7) 



So in this case, 7 = ±r lies on the boundary of the accessible region, 7 G 
dUh,i- In case (i), 7 lies inside, in case (ii) outside. This will be relevant for 
the topology of V\ /5 see next section. 

Let us now consider 7 which are not collinear with r. We show in the 
Appendix that given 7, the preimages in V\ l are the points (7, 1) with 



(7,^7) V ( v , A ~ lv ) 

where v = ( J 4 _1 (r x 7)) x 7. This tells us that for all 7 with Ui(j) < h, or 
7 from the interior of Uhj, there are exactly two preimages in V\ v For 7 
on the boundary, 7 G c^P^j, the preimage is unique. No real J exists for 7 
outside Uh,i- 

Computing the determinant A (7, 1) according to Eq. (11), we find 

A( 7 ,Z) = (A- 1 ix 7 ,^ 1 (7xr)). (15) 

This is zero for 7 G dUh,i, where A~ x l x 7 = 0, corroborating the assertions 
made in connection with the right panel of Fig. 2: the parts of V\ x where 
A > and A < each project to the entire accessible region U^x C S 2 (7). 
(Note that Uh,i , hence also E\ x and V\ l , may consist of several disconnected 
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components.) The other possibility for A to become zero, 7 x r = 0, is of 
another kind; its projection to S 2 (7) gives at most the two points 7 = ±r. 
If they are inside Uh,i, their preimages in Vi, are circles, see case (i) above 
Eq. (13). The consequences are discussed in the next section. 

The condition S = for tangency of trajectories with V\ t is obtained 
from (12) by differentiation: 

S = (A~ l (l x A~H - 7 x r), r x 7) + (A'H, r x (7 x A~H)} = . (16) 

This must be considered together with (14). We do not attempt to derive an 
explicit expression for the projection of these lines to the Poisson sphere, but 
they will be shown in the graphical representations of the following sections. 

4 Topology of V\i and its projections 




Figure 4: Projection of V\x onto the Poisson spheres S (7) (left) and S + (7) 
(right). On the abscissa tp goes from to 2ir, the ordinate •& varies from (top) to 
7T (bottom). A= (2,2,1), r= (1,0,0), ft = 3.4, Z = 2.87. 

Fig. 4 shows a typical projection of V\ 1 to two copies S~ (7) (left) and 
S + (7) (right) of the Poisson sphere, for A = (2,2,1) and r(l,0,0). The 
coordinates of the two panels are < ip < 2ir on the abscissa and < § < tt 
on the ordinate. The white regions at small and large 1? are inaccessible with 
the given values of ft = 3.4 and I = 2.87. The grey and black parts together 
are the accessible region Uhi which has the topology of an annulus. In the 
black parts the trajectories are incoming, S < 0, in the grey region they are 
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outgoing, S > 0; the boundary between these regions is the set of tangency. 
If we were to record only incoming trajectories, only the black parts needed 
to be considered, but note it appears on both copies S =F (7) . 




Figure 5: For a topologically more adequate representation of the preimages in 
V\ i of the poles 7 = ±f , these poles are replaced by "polar circles" which turns 
the Poisson sphere into a cylinder. 

The white dots at {&,ip) = (ir/2,n/2) are the "north pole" 7 = r, the 
white dots at (tt/2,3tt/2) are the "south pole" 7 = — f. We know from the 
previous section that these points are special. When they are physically 
accessible as they are here, then their preimage in V\ j is a circle (instead of 
a point as for all other 7). This suggests to represent V\i not in projection 
to a sphere but rather to a cylinder: punctuate S 2 (7) at the poles and insert 
"polar circles" there, as illustrated schematically in Fig. 5. Instead of the 
two spheres we then have two cylinders for the projections of and 
V^i- But notice that these two cylinders share the same circles at the poles, 
hence together they form a torus T 2 which we shall call the "PP-torus" 
T 2 (7), the two P's referring to Poisson and Poincare, or to the two copies of 
the punctuated Poisson sphere, and the subscript 2 indicating that each 7 
is represented twice. We denote its two halves by T~, T + , and remark 
that each carries one copy oiU^i. and V^i, respectively, project 1:1 to 
interior points of the two Uh,i, and their boundary, defined by A(7,Z) = 0, 
projects to dUhj- Thus, if we identify the two copies of dUh,i on T~ and 
T + , we obtain a 1:1 map Vii — * T 2 (7) which preserves the topology. 

The graphical implementation of this idea is automatically achieved if 
we first transform to new polar coordinates (#', ip') where 7 = r is taken as 
polar axis instead of 7 = (0, 0, 1), and then interpret = and $ = 7r not 
as points on a sphere but as circles on a cylinder. An example is shown in 
Fig. 6 where the two representations are compared. The top row exhibits 

as in Fig. 4, the bottom row T^. The north pole = (7r/2,7r/2) is 

transformed into the circle = (upper boundary) , and the south pole into 
the circle = it (lower boundary). Identification of these circles between T~ 
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and T + produces the full torus on which the non-accessible white regions 
are now the two holes in the neighborhoods of the points (■#', ip') = (vr/2, 0) 
and (7r/2,7r). Since corresponding boundaries dUh,i on T~ and T + must 
also be identified, each pair of holes on the two sides generates a handle. 
Hence, the surface of section V\ l for this example has the topology of a 
torus with two handles, i.e., it is a manifold M|. The lines of tangency, in 
this particular example, seem to have become lines of constant . 




Figure 6: Comparison of the projection of V\ z onto two copies S T of the Poisson 
sphere S 2 (#, ip) (top) and onto the two halves T T of the torus T 2 ^??',?//) (bottom). 
The parameters are A = (2,1.5,1), r = (1,0,0), h = 80.5, I = 12.80. Incoming- 
intersections of the pink orbit are all on the left, outgoing on the right. 

One trajectory has been added, in pink color. Its incoming part S < 
lies entirely on T~, the outgoing part S > on T + . Note that the incoming 
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part of the trajectory comes close to the pole 7 = r whereas the outgoing 
part has points close to the south pole 7 = — f . The comparison of upper 
and lower row illustrates that the torus T| is the most natural basis for a 
projection of V\ ; . 

Simpler manifolds occur, e.g., when Uh,i is a topological disk on S 2 (7). 
If this disk does not include any of the poles ±r, its representation on the 
PP-torus is two disks, one on each half. The surface of section is then the 
union of two disks which are glued together at their boundary; in other 
words, V\ { is a topological sphere S 2 . If the disk contains just one of the 
poles, its representation on the torus is a cylinder, and V\ l is the torus T 2 
obtained by identifying its upper and lower rim. If both poles are contained 
in the disk, the same kind of arguments show that V\ l is a 2-manifold of 
genus 2. 





Topology 


f a connected component of £^ t 


MP 3 


s 3 


S^S 2 


(S 1 xS 2 )#(S 1 xS 2 ) 


Structure of Uh,i 


| number of poles in Uh,i 


S 2 


D 2 


S 2 \2D 2 


S 2 \3D 2 







S 2 




M 2 


1 






M 2 


M 2 


2 




M 2 


M 2 


M 2 



Table 1 : Possible topologies of the Poincare surface of section V\ 



It is known from the work of Tatarinov [11] and Bolsinov et al. [2] that 
connected components of Uhj C S 2 (7) come only in four topological kinds: 
the entire sphere S 2 , a disk D 2 = S 2 \D 2 , an annulus S 2 \2D 2 , and a sphere 
with three holes S 2 \3D 2 . The corresponding topological types of V\, are 
listed in Tab. 1, for all possible locations of the poles relative to Uh,i- They 
are manifolds of genus g with g ranging from to 4. If we distinguish pairs 
of £\ l and V\ l when either of the partners is different, then the table shows 
there are 10 different possibilities. 

To wrap it up, we propose to start with the PP-torus T^??',^') as the 
appropriate manifold onto which the Poincare surface of section V\ l is to be 
projected 1:1. It is constructed from two copies of the Poisson sphere where 
+r serves as reference point for the polar coordinates (i?', ip'). The spheres 
are punctuated at the two poles +f and — r which are replaced by circles 
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(#',^/) = (0, if}') and (it, ip'), respectively. Identifying corresponding circles 
on these two cylinders produces the torus. It will be seen in the next section 
that this identification involves a relative shift of the two angles ip' by ir. 
The PP-torus carries two copies of the accessible region Uh,i- Boundaries 
dUhj of these two copies are to be identified. The resulting manifold is a 1:1 
image of the complete Poincare surface of section V\ l . 



5 Examples 

In this last section we give nontrivial examples of surfaces of section V\ l 
where the extra bifurcation lines (13), or 

h - (17) 

do not coincide with bifurcations of the energy surface. This cannot happen 
in the Katok family of systems considered in [6] (which includes all cases 
of Lagrange and Kovalevskaya) because there the center of mass r lies on 
a principal axis, and then (17) is the bifurcation line corresponding to the 
relative equilibrium where the body rotates about that axis. In order to 
make (17) different from any of the relative equilibrium conditions, we must 
choose r off the principle axes. An arbitrary example is shown in Fig. 7, 
with A = (2,1.1,1) and r = (0.94868,0,0.61623). The two orange lines 
separating from the red and blue relative equilibria correspond to 7 = r 
(left) and 7 = — f (right). When h is increased at fixed I, Uhj is at first a 
disk which does not contain any of the poles. Upon crossing the left orange 
line, +r enters Uh,u then at the right orange line, — r comes in as well. 

The six lines in Fig. 7 show a partition of the (h, £)-plane (with I > 0) 
into eight regions (plus the non-physical region to the left of the red line). 
Seven of them are represented by the black dots on the line I = 3.25. The 
tiny eighth region lies between the yellow and green lines to the left of the 
orange line (black dot at (h, I) = (3.1486,2.72)). An additional region, not 
contained in the figure, exists because at large values of h and I the left 
branches of the yellow and orange lines intersect; point (h, I) = (14.0, 7.2) is 
representative for this region where the left orange line lies to the right of 
the yellow line. 

These nine regions correspond to nine different topologies of V\ l: includ- 
ing six of the types listed in Tab. 1, and three examples of disjoint surfaces 
of section. We discuss them one by one. The following Poincare sections 
contain typical orbits, with different colors for different orbits. When certain 
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orbits are the same type (elliptic or hyperbolic) in different regions they are 
given the same color. 



1 




3.142 3.149 3.156 

h 



Figure 7: Bifurcation diagram for the Poincare surface of section V\ t , for A = 
(2,1.1,1) and r = (0.94868,0,0.61623). The red, blue, green, and yellow lines 
correspond to critical values of (7), the two orange lines come from (17). Top right: 
Magnification of the rectangle marked on the left. Bottom left: Magnification 
around the point (h,l) = (3.1486,2.72). 

The first example (h,l) = (1.8,3.25) represents the leftmost region be- 
tween the red and orange bifurcation lines. It is shown in the top row of 
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Fig. 8. The accessible region Uh,i on T T is a rather small topological disk 
which does not contain the center of mass r. Gluing the two copies together 
along the boundary oiUh^i we obtain for V\i a topological sphere S 2 . 

The second row of Fig. 8 shows the case (h, I) = (2.6,3.25). The north 
pole (upper rim) has entered Uh,i- The orbit structure reveals that when 
we identify the two circles corresponding to f on T T , the angle ip' must be 
shifted by it. On the Poisson sphere, Uh : i is again a disk D 2 , but on the 
torus T + UT~ its two copies together form an annulus. Gluing together the 
two boundaries dUh,i on T 2 , we find that V\i is a torus T 2 . 

Increasing the energy beyond the yellow line in Fig. 7, the topology of 
£? , changes from S 3 to 2S 3 , a new disjoint disk appears as part oiUhi- This 
disk contains no pole, hence we deduce that V\ , is the union of a torus 
T 2 and a sphere S 2 , see the third row of Fig. 8 with (h,l) = (3.8,3.25). It 
may be surprising that the trajectories shown indicate only regular behavior 
with no sign of chaos even though the equations of motion are certainly non- 
integrable. 
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Figure 8: Projections of V\ l onto T (left) and T + (right). Parameters are those 
of Fig. 7. The angular momentum is I = 3.25 for all cases, the energy from top to 
bottom is 1.8, 2.6, 3.8 and 4.45. 



18 



The next point (h,l) = (4.45,3.25), see the right part of Fig. 7 and 
the last row of Fig. 8, represents a region where both disjoint disks of Uhi 
contain a pole. The surface of section is therefore 2T 2 . It appears that the 
accessibility of the south pole (corresponding to the center of mass in upright 
position) suddenly introduces manifest chaos. While the component around 
r still looks very regular, the one containing — r is dominated by chaotic 
motion. The orange and purple orbits belong to the same chaotic region but 
for a long time remain separate, indicating the presence of cantori related 
to the chain of islands between them. 

The next three cases with / = 3.25 are shown in Fig. 9. The first row 
with h = 4.675 represents the narrow (h, Z)-region between the green and 
yellow bifurcation lines, see the right part of Fig. 7. Uh,i is the Poisson sphere 
minus two disks, i.e., a topological annulus; as this annulus contains both 
poles, the surface of section V\i is a manifold M| of genus 3. Compared to 
the previous case h = 4.45, the two disks about the poles have merged at two 
points and now allow the chaos to sneak into the upper region. Nevertheless, 
the neighborhood of f remains mostly regular. 

One of the two inaccessible disks disappears as h is increased beyond the 
yellow line in Fig. 7, see the second row of Fig. 9 where h = 5.3. Uhi 1S now 
a disk which contains both poles, hence V\i is a manifold M| of genus 2. It 
appears that the phase space is fairly distinctly divided into one chaotic and 
two regular parts. The "outgoing" intersections of the chaotic orbits take 
place in a region connected to the south pole r, the "incoming" lie between 
the poles. A bigger regular region surrounds the north pole, a smaller one 
connects to the south pole at about one half of the ^'-circle. So chaos and 
order meet at the south pole, depending on the angle ip'. 

Finally, when h crosses the blue bifurcation line, the whole Poisson 
sphere becomes accessible. This is shown, with h = 6.5, in the last row 
of Fig. 9. The surface of section is now isomorphic to the entire PP-torus 
T|. As to the distribution of regular and chaotic orbits on this torus, we 
observe again two regular parts separated by two chaotic regions. 
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Figure 9: Projections like in Fig. 8. Again I = 3.25 in all three cases; the energy 
from top to bottom is 4.675, 5.3 and 6.5. 
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Two more types of surfaces of section are shown in Fig. 10. They do not 
exist for I = 3.25. The upper row shows the example (h,l) = (3.1486,2.72); 
the corresponding dot in Fig. 7 lies between the green and yellow lines, but 
to the left of the orange line. Compared to the first row of Fig. 9, the lower 
part does not contain the south pole, yet there remains a small hole (barely 
visible in the figure) inside the accessible region, near the center. Uh,i is 
a topological annulus containing the north pole, hence V\i is of type M|. 
In the scenario where this case emerges from the third row of Fig. 8, an 
outbreak of chaos takes place in the lower part. 

The last case is (h,l) = (14.0,7.2); the accessible region consists of two 
disks neither of which contains a pole. Hence V? r is the union of two topo- 
logical spheres. 



21 



With this series of nine types of Uh,i we have illustrated 6 out of the 10 
possible topologies listed in Tab. 1. But notice that our choice A = (2, 1.1, 1) 
represents only one of the seven different classes of moments of inertia that 
Katok distinguishes with r = (1,0,0) (the class K2 in the notation of [6]). 
With other choices it is possible to also find examples for the other types of 
U Kl . 

6 Summary and outlook 

We have presented a new tool for the study of rigid body dynamics: a com- 
plete Poincare surface of section V\ { and its 1:1 representation in projection 
to a suitably defined torus, the PP-torus T|. The surface V\ x C £f ll is 
defined by local extrema of the vertical component of the center of mass r, 
d(7,r)/dt = 0, and the PP-torus is obtained from two copies S + (7) and 
S~(7) of the Poisson sphere by the following procedure: punctuate S ± (7) 
at the points 7 = ±r/r, replace these points by "polar circles" to obtain two 
cylinders T =t (7), and then identify corresponding polar circles on T + (7) and 
T~ (7) . The torus T| so constructed contains two copies of the accessible 
7-region, and V\i projects 1:1 to the union of these two copies. Identifying 
the boundaries dUh,i on T + (7) and T _ (7) produces a manifold which is 
homeomorphic to V\ , and at the same time readily accessible to intuition. 

We have shown examples for nine topologically different surfaces V\ 1 
and propose to use the tool for further exploration of the dynamics of rigid 
bodies in their 4-dimensional parameter space. In a forthcoming publication 
we will apply it to the family of systems A = (2, 2rj, 1), r = (1, 0, 0), which 
contains the integrable system of Lagrange (77 = 1) and Kovalevskaya (77 = 
2). Pictures like those of Figs. 8-10 reveal the fate of particular features of 
the phase space structure under parameter variation, such as location and 
stability of isolated periodic orbits, their bifurcation schemes, as well as the 
extent and entanglement of regular and chaotic regions. 

A Appendix: Explicit characterization oiV\ l 

Given an energy h, an angular momentum I, and a point 7 on the Pois- 
son sphere (not collinear with r), we determine the corresponding vectors I 
on the Poincare surface defined by S = 0, or (i,^4 _1 (r x 7)) = 0. For 
abbreviation we introduce the notation 

u := A _1 (r x 7), v := u x 7 . (A.l) 
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The equations S = (l,u) =0 and L = (l,~f) = I define planes in /-space 
which intersect in a line with direction v. It is easy to check that 

(v,v) 

is one point on this line, hence we may parameterize it with I = lo + vt. The 
energy equation (3) determines the possible values of t. We write it in the 
form 

(lo + vt, A' 1 (l + vt)) = 2(h + (r, 7 )) =: c. (A.3) 
This gives a quadratic equation for t, 

(v, A^v)^ + 2(l , A^vjt + <i , A- 1 ^) - c = . (A.4) 

Inserting its solutions into I = Iq + vt, we find I = l\ + I2 with 

(loiA^v) 

h=lo--, T3T v v ( A - 5 ) 

(v,A i v) 

and 

= {v J-i v) ^J(h,A-±v)i - ((l ,A-Ho) - c)(v,A-±v). (A.6) 

To evaluate further, we use an identity which holds for arbitrary symmetric 
matrices A and vectors a, b: 

AaxAb = (detA) A' l (a x b). (A.7) 

For example, with v = u x 7 we find = (Au x J 47)/detA and 

(«, A~ l v) = — xj,Aux Aj) = _L( Uj Au)(7,A7) ; (A.8) 

in the last step we used Lagrange's identity and (7, Au) = (7, r x 7) = 0. 
With similar arguments, and using the explicit form (A.2) for Iq, we find 

h = 1 1 ~~t r • (A.9) 
(7, A-y) 

To evaluate I2, we first employ Lagrange's identity to obtain 

(l Q ,A~ l v) 2 - (I^A'Hq)^^-^) = (IoXv^^vxA'Hq), (A.10) 
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and then (A. 7) with A 1 instead of A to transform this into 



-1 



(l x v,A(l x v)} . 



(A.ll) 



detyl 



An easy calculation using (v,u) =0 shows that Iq x v = lu. Inserting this 
into the radicand of (A. 6) we find that it may be written as 
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